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^ ■ ABSTRACT 

> : 

T-H ■ We present results from simulations on two aspects of quark confinement in the 

I pure gauge sector. First is the calculation of the profile of the flux tube connecting 

a static qq pair in SU (2). By using the Michael sum rules as a constraint we give 
evidence that the energy density at the center of the flux tube goes to a constant 

■ as a function of quark separation. Slow variation of the width and energy density 
\ is not ruled out. Secondly in the confined phase of lattice U(\) we calculate the 

■ curl of the magnetic monopole current and show that the dual London equation is 
I ' satisfied and that the electric fluxoid is quantized. 

D . 
■ 

• • ■ 1. Introduction 

> ■ 

•i-H . 

^ ' I would hke to report on efforts of the Cracow-LSU collaboration to study the 

■ response of gauge fields to the presence of static sources. Confinement is the central 
issue which is seen as a consequence of the vacuum squeezing the field lines to form 
flux tubes connecting color charges. One is beginning to see considerable detail of 
the field distributions in the flux tube, e.g. its size and shape, the chromoelectric 
and chromomagnetic field components, and monopole currents that are responsible 
for squeezing the tube to form an Abrikosov vortex. 

In the first part of this talk I will describe our results for SU(2) lattice gauge 
theory [1, 2, 3]. We examine the flux tube between a qq pair, calculating the six 
chromoelectric and chromomagnetic components to the energy density and action 
density. We parametrize the profile of the flux tube and study scaling. The results 
are subjected to the check provided by the Michael sum rules[4]. 

The second part of this talk concerns the mechanism that leads to flux tube 
formation. This is the role of the solenoidal magnetic monopole currents that surround 
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the flux tube. In this work[5] we show in the confined phase of U(l) that the curl of the 

monopole current has a profile similar to the electric field and that the dual London 
equation is satisfied and electric fiuxoid quantization occurs. We demonstrate that the 
the fiux tube is precisely the dual of the Abrikosov vortex in Type II superconducting 
materials. 



2. Flux Tubes in SU(2) ^ 

2.1. Background 

In this calculation we measure the field energy densities by correlating the small 
plaquette with the Wilson loop. Full details of the simulation are given in Ref. [2], 
which also contains further references. Further details of the fiux profiles will appear 
in a companion paper Ref. [3]. By fixing our attention on the middle time slice of 
the Wilson loop, the time-like segments form world lines that approximate a static qq 
pair. The 3 space-space plaquettes measure the magnetic component of the energy 
density and similarly the 3 space-time plaquettes measure the electric components. 

Before defining the flux calculation in more detail, we point out that the Wilson 
loops themselves are used to extract the transfer matrix eigenvalues which give the 
static quark potential and are further used to extrapolate the flux measurements to 
infinite time extent of the Wilson loops. Specifically we determine the eigenvalues of 
the transfer matrix by fitting the Wilson loops to the exponentials as described in 
Ref. [2]. 

< W{R, T) >= ^ Ae-''^^^)^. (1) 

i 

Eq{R) is of special interest since it contains the static quark potential: 

E.(H) = -^ + .«+£||. (2) 

The term independent of R is the self energy of the two quarks which does not scale 
but diverges as a ^ 0. Since Wilson loops can be calculated quite accurately, the 
static potential is a useful physical quantity to check scaling and thereby determine 
the lattice spacing a{P). All our data is consistent with standard values a(2.3) = 0.171 
fm, a(2.4) = 0.128 fm and a(2.5) = 0.089 fm. 

2.2. Flux Tube Profiles 

The lattice observable needed to measure the flux is the following [6, 1, 2]. 
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where W is the Wilson loop, P^'^ the plaquette located ai x, P — and xr is a 
distant reference point. In the classical continuum limit 

r ^-^ -l{iFn\,-vac, (4) 

where the notation (■ • ■)qq-vac means the difference of the average values in the qq 
and vacuum state. Prom now on we shall be using field components in Minkowski 
space and hence 

r - li-Bl -Bl -Bl El El El). (5) 

Correspondence between various components and f^" is standard: space-space pla- 
quettes are magnetic, space-time plaquettes are electric. The energy and action den- 
sities are respectively 

e=\{E^ + B% 

^^\{E''-B% (6) 

Since the magnetic contribution turns out to be negative, there is a strong cancellation 
between the two terms in the energy, whereas they are enhanced in the action. 

Figure 1 gives the flux profiles. The cancellation which suppresses the energy 
density in the flux tube is evident. However notice that the self energy of the quarks 
is not similarly suppressed. This follows because the self energy is primarily electric. 
We fitted the energy and action density in the plane at the midpoint between q and 
q using the function 

/(r^) = ae^v{-\lh^ + {rjcf). (7) 

The peak value and the width at half maximum were very well determined using a 
fit for each of 70 cases of different loop sizes and values of (3. For the third parameter 
we chose the decay length of the tail of this function and found it less well determined 
but with a value typically close to the width at half maximum. The details of the 
analysis will be given in a forthcoming paper [3]. Here we just give the results of the 
extrapolation to infinite Wilson loop time extent in Figs. 2 and 3. 

The basic issue is whether the peak value of the energy density stabilizes to a 
constant or goes to zero with quark separation. This is not easy to settle as can be 
seen in Fig. 2. Roughly speaking we know from the linearly rising potential that the 
string tension ~ (width)^ x (peak value) should be constant. Both Figs. 2 and 3 
show that we are marginally asymptotic in quark separation, R. The two curves are 
~ and ~ A Coulomb field would fall hke the former but since the string 
tension is constant the asymptotic width would have to grow like which clearly 
it docs not. Therefore we can rule out a Coulomb field as expected. Interestingly 
for small separations, the eyeball fit to ~ is quite good which may be due to a 



Figure 1: Energy and action profiles: (a) energy density in the plane containing qq; 
(b) energy density in the plane midway between q and q; (c) and (d) similarly for 
action density. 



Figure 2: (a) Peak value of energy and action density; solid squares: p = 2.5; triangles: 
P — 2.4; open squares: /3 = 2.3. The two curves are 1/R and 1/R'^ arbitrarily 
normalized; (b) blowup of (a). 



Figure 3: Width at half maximum for energy and action density. 



Coulomb like behavior at small distances. (The above argument that the width must 
grow like does not apply because there is no string for small R.) The dielectric 
model[16] predicts the peak density ~ 1/R. This function (arbitrarily normalized) 
does not seem to fit the data very well. However such a behavior would imply the 
width ~ \/R which is certainly possible in OTir data. We can say that the peak energy 
density and width are consistent with a constant value for large quark separation but 
we can not rule out a slow variation. The issue can be tightened by making use of 
the Michael sum rules[4] as we mention in the next section. 

Figure 4 illustrates a general feature of our data. The cluster of three points for 
each R and T correspond to the three quantities: 

e - liEl + BD + liEl + Bl), 

e{T^R) = l{El + Bl)-\{El + Bl). (8) 

If one turns the Wilson loop on its side, the || components are unchanged but the 
_L components of the electric and magnetic fields are reversed: E\ <-> —B\. Hence 
there is a sign change in the third expression. The central points of the cluster are 
the II components only. The clustering of the points implies that the ± components 
of the electric and magnetic contributions to energy density are approximately equal 
hut of opposite sign and cancel. The width of the peak is even less sensitive to the 



Figure 4: Peak density for RxT Wilson loop sizes R — 3 — 6, T = 3 — 6. For fixed T 
the points decrease monotonically with R. Triangles: energy density for R x T loop; 
circles: energy density for T x i? loop; squares: || components of E'^ and only 

transverse components giving essentially the same value for all three points. 
2.3. Sum Rules 

A consistency check on the flux distributions can be obtained by using the Michael 
sum rules[4] for energy and action. 

l^iEixf + Bixf) = Eo{R), 

X 

WiE{xr-B{xr) = -p^[E,{R)-'M]_pm. (9) 

2 ^ a a a 

X 

Here Eq[R) is given by Eqn.(l), and (■ = ^). In rcf. [2] we have shown that our data 
are essentially consistent with these sum rules. The one difficulty is the fact that 
the self energy, c(/?)/a(/3), determined from the potential differs from the self energy 
determined from the the action sum rule. This may be due to an ambiguity in the 
definition of self energy or possibly due to our classical expressions for energy and 
action which ignores quantum corrections. By taking a derivative of these expressions 
with respect to the quark separation, R, this difficulty is avoided. This gives the 
relation 

<7^ = -A; aA^\Y.iE{xf-B{Sn a ^ ^ + ). (10) 

The sums are now over the plane midway between the qq pair. 



Using the sum rules we find that the f3 function {—(5d/a) ^ 10. ±2. compared to the 
current estimates 7.±1.. The asymptotic value is — 51/121 + 37r^/5/ll = 6.0(/? = 2.4). 
There is ample evidence from other measurements that although scaling works well, 
asymptotic scaling is violated[8] and hence we do not expect to get the asymptotic 
value. 

An alternative approach is to assume the sum rules are correct and use them to 
infer information about the energy density from the action density which is far easier 
to measure since relative errors are down by an order of magnitude. As is clear from 
the sum rule, the action does not scale yet the variation over these values of P is 
very small. An examination of Fig. 2(a) shows that the action for each /3 seems to 
stabilize to a constant for increasing distance for the peak density and for the width. 
This is quite striking for f3 = 2.3 and 2.4. For = 2.5 R appears to be too small 
to draw a conclusion. These data do not suggest that the peak value is tending to 
zero at all. We would like to use the sum rules to predict the behavior of the energy 
density. A constant peak energy density follows only if the widths of the energy and 
action peaks have the same behavior. Figure 3 shows that in fact they do. From this 
and using the sum rules we conclude that the energy density stabilizes to a constant 
value also. This conclusion is an argument against the dielectric model[16]. However 
we have little to say about logarithmic behavior of the fiux tube width as predicted 
by Liischer[7]. For more details see Ref.[3]. 

3. Mechanism for confinement in U(l)'^ 

3.1. Dual Superconductor Model of Confinement 

We now turn to the mechanism for flux tube formation and present direct evidence 
that supercurrents of magnetic monopoles produce a dual Abrikosov vortex[9]. U(l) 
lattice gauge theory in 4 dimensions has both a confined phase at large charge and a 
weak coupling deconfined phase corresponding to continuum electrodynamics with a 
Coulomb interaction between static charges. Therefore confinement or its absence can 
be studied using U(l) lattice gauge theory as a prototype, before tackling the more 
complicated non-Abelian theories that actually describe quarks. Much evidence for 
the dual superconductor hypothesis has accumulated from studies[10, 11, 12] of lat- 
tice gauge theory. PolyakovflO] and Banks, Myerson and Kogut[ll] showed that U(l) 
lattice gauge theory in the presence of a quark-antiquark pair could be approximately 
transformed into a model describing magnetic current loops (the monopoles) inter- 
acting with the electric current generated by the qq pair. DeGrand and Toussaint[12] 
demonstrated via a numerical simulation that the vacuum of U(l) lattice gauge the- 
ory was populated by monopole currents, copious in the confined phase and rare in 
the deconfined phase. This behavior has also been seen in non-Abelian models after 
gauge fixing[13]. Many studies of non-Abelian models using Dirac monopoles[13, 14] 
or other topological excitations [15] support the dual superconductor mechanism, al- 
though other studies[17] dissent. 
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So far, studies of confinement have examined "bulk" properties such as the monopole 
density[12, 13], and the behavior of the static quark potential[14]. In a recent paper[5] 
we presented the first direct evidence that the flux tube is a dual Abrikosov vortex. 
We further show that there are exact U(l) lattice gauge theory analogues of two key 
relations that lead to the Meissner effect in a superconductor; the London equation 
and the fluxoid quantization condition. 

3.2. Electric Field Profiles 

Our simulations were done on a Euchdean spacetime lattice of volume 9^ x 10. 
The static charges are represented by a Wilson loop as in the previous section. We 

take a 3 X 3 loop in the z — t plane and measure the fields in the x — y plane 
at the midpoint between the charges. Because of the geometrical symmetry of the 
measurements only the 2;-components of {£) and (V x Jm) are nonzero. If the Wilson 
loop is removed, even the z-components average to zero, so the response is clearly 
induced by the presence of the static charges. Only the imaginary part of the Wilson 
loop contributes to the averages of these two quantities. 

The plaquette measures flux passing through a unit square on the lattice^ 

exp[iea2F^,(r)] = exp[z^^,(r)] = U^{f)U,{f+ fi)Ul{r + iy)Ul{f). (11) 

The electric flux in lattice variables is 

S^{f)=lmexp[ie^,{f)]. (12) 

Figure 5(a) shows the electric flux distribution for /? = 1.1 where the vacuum is 
in the deconflned phase. The broad flux distribution seen is identical to the dipole 
fleld produced by placing two classical charges at the quark positions, except that the 
classical value of the flux on the qq axis is a factor of two smaller. We measure the 
total electric flux from one quark to the other, including not only the flux through the 
plane between the charges (0.8504 ± 0.0045) but also the flux (0.0951 ± 0.0028) that 
flows through the lattice boundary because of the periodic boundary conditions. This 
yields a total flux of (0.9453 ± 0.0053), close to the theoretical value $e — e/Vhc — 
1/^ = 0.9534. 

Figure 5(b) and 6(a) shows the electric flux in the conflned phase {f3 = 0.95). In 
this case the flux is conflned almost entirely within one lattice spacing of the axis and 
essentially no flux passes the long way around through the lattice boundary. The net 
flux is again equal to l/-\/^ within statistical error. This behavior is exactly what one 
would expect from the superconducting analogy, where the flux has been "squeezed" 
into a narrow tube. 

3.3. Magnetic Monopole Supercurrents 

The monopole currents are found by a prescription devised by DeGrand and 
Toussaint[12], which employs a lattice version of Gauss' Law to locate the Dirac 

^The flux here means Ex (area). We use the same term for E"^ or since it has become an 
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Figure 5: Surface plot of the electric flux through the xy plane midway between the 
qq pair when the system is in (a) the deconfined phase {(3 = 1.1) and (b) the confined 
phase {P — 0.95). The hne joining the pair is located at (0,0). 

string attached to the monopole. The net flux into each plaquette face is given by 
(^^zy(^) mod 27r). If the sum of the fluxes into the faces of a 3- volume at flxed time 
is nonzero, a monopole is located in the box. A non-zero net flux can occur only if a 
multiple of 27r arises from the mod operation. Therefore the net flux is 

a^Fni, = n{2TThc/e) = ucm (13) 

Gfaces 

The net flux into the box at flxed time thus yields the monopole "charge" density, 
the time component of the monopole 4-current Jm- The spatial components are 
found similarly. The monopole currents form closed loops due to the conservation 
of magnetic charge. Finally the curl is calculated by the line integral of the current 
around a dual plaquette. 

— * — * — * 

We show in Fig. 6(a) (S) and in Fig. 6(b) — (V x Jm) in the conflned phase as a 
function of the distance from the qq axis. The data show that the spatial variation of 
the flux and the curl of the current are very similar, except for the point on the axis 
which will be discussed below. 

3.4- London Equations, Fluxoid Quantization and the Abrikosov Vortex 

In order to interpret this result we would like to review the ordinary London 
theory. In the next section we interchange electric and magnetic quantities to get the 
dual results. 



Figure 6: Behavior of (a) the electric flux, (b) the curl of the monopole current and 
(c) their sum, the fluxoid, in the confined phase (/? = 0.95) as a function of the 
perpendicular distance r±_ from the qq axis. The solid line in (a) gives the flux for a 
pure Coulomb field . 



A concise statement of the London theory is contained in the relation^ 



A+—J^O; (V-A = 0). (14) 
c 

If the charge density is zero, then in this gauge the electric field is given by —A/c and 
therefore E = \^J/c. This describes a perfect conductor and is just Newton's law for 
free carriers, eE = mv. By taking the curl of Eqn.(14) we obtain the condition for a 
perfect diamagnet. 

V X J = --^B. (15) 

This relation together with Ampere's law V x B = J/c gives V^-B = B/X^ which 
implies that the magnetic field falls off in the interior of the superconductor with a 
skin depth A. 

Finally the fluxoid is given by the integral 

/ (B + —V xJ)-nda^ / (A + —J) ■ dl^ n— = (16) 
Js c Jc c 2e 

If the curve C is in a simply connected region of a superconductor, then n = 0. 
However if the curve encircles a hole in the material then n need not be zero but must 
be an integer. In an extreme type II Abrikosov vortex, a very small core is comprised 
of normal material. A single unit of magnetic flux of radius ~ A passes through the 
vortex. The fluxoid density 5 + x J is zero everywhere except in the region of 
the normal material. In the limit in which the core is a delta function we obtain: 

A^ 

B + —V X J^^rnS^(x±)n,. (17) 
Further if we use Ampere's law we can get an analytic expression for B^. 

B,{r^) - X^V^B, = ^m8\x^)- B, = ^K^{rJ\). (18) 



3.5. Dual Superconductor 

We interpret our results using the following relations which are the dual of the 
corrsponding relations in the previous section: 

E-^VxJr^^ $e5'(Xi)n,. (19) 

where the unit of electric flux $e = l/V/^- Figure 5(c) shows the result of fltting our 
data to this relation. The London penetration depth. A, is the only free parameter. 
We are able to determine a value of A which gives zero away from the axis and a delta 

^We use Heaviside-Lorentz units to be consistent with lattice gauge theory. 



function with the correct coefficient on axis. Further we can check that the electric 
flux profile is given by 

E. = ^^o(r/A) (20) 

This function has no free parameters. This curve is also shown in Fig. 6(a) showing 
excellent agreement. Wc find a value of A/a = 0.482 ± 0.008, which is consistent 
with the range of penetration of the electric flux in Fig. 6(a) and the thickness of the 
current sheet in Fig. 6(b). We expect that, as in a superconductor, the transition 
to the deconflned phase will be signalled by a divergence of the London penetration 
depth. We have therefore measured A further from the deconflnement transition at 
13 = 0.90, and flnd a smaller penetration depth of A/a = 0.32±0.02. In the deconflned 
phase, j3 = 1.1 wc flnd an almost insigniflcant value of (V x Jm) and fltted values of 
A were larger than our lattice size. 

In summary, one can flnd a value of the London penetration depth. A, that satisfles 
Eqn.(19) off axis. One then flnds that for the point on axis, the same value of A gives 
one quantum of electric flux as predicted by Eqn.(19). Finally the same value of A 
gives a good flt to the proflle using Eqn.(20). Hence considerable detail of the dual 
Abrikosov is verifled. It is perhaps surprising that a nonlinear, strongly interacting, 
model such as U(l) lattice gauge theory could be described by such a simple model 
as the linear London equations but our results indicate that the operators {£) and 
(V X Jm) when measured in the presence of source of external flux like a Wilson loop, 
give an unambiguous indication of the conflnement of electric flux by a monopolc 
current distribution. The simulation yields a large signal even with modesnt amounts 
of computer time on a Sun workstation. Although the Meissner effect itself requires 
only that Eq. (20) hold off axis, our data also support the more restrictive fluxoid 
quantization relation on axis. This additional relation reflects the single- valued nature 
of the order parameter in a Ginzburg-Landau description of the monopole condensate. 
Because the monopoles appear pointlike in our simulations, lattice gauge theory looks 
like an extreme type-II superconductor. 
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